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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 


An  approximate  solution  is  given  for  the  effect  of  an  exponentially  decaying 
pressure  pulse  traveling  with  superseismic  velocity  on  the  surface  of  a  half 
space.  The  material  of  the  half  space  is  an  elastic-plastic  model  of  a 
material  having  internal  Coulomb  friction.  The  yield  condition  selected  may  be 
suitable  for  a  granular  material.  The  effect  of  a  step  wave  for  this  geometry 
and  medium  was  treated  previously.  For  that  case,  the  peak  pressures  do  not 
decrease  with  increase  in  depth,  while  such  a  decrease  is  obtained  for  a 
decaying  surface  load.  It  was  the  prime  purpose  of  this  investigation  to 
determine  the  magnitude  of  this  attenuation.  The  approximate  solutions  obtained 
are  valid  for  a  limited  distance  behind  the  wave  front,  and  are  tabulated  for 
11  different  sets  of  parameters  pertaining  to  the  material  and  velocity.  The 
tabulated  results  show  that  the  peak  pressures  in  the  case  of  the  decaying 
surface  load  do  decrease  with  depth,  but  that  the  decrease  is  less  than  one 
might  Intuitively  expect. 
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SECTION  I 


INTRODUCTION 


In  a  preceding  report,  Ref.  [l]*^,  the  effect  of  a  step  pressure,  Fig.  1, 


progressing  with  superset smic  velocity  V  on  the  surface  of  a  half -space  has 


been  studied  for  the  limiting  case  k  •*  0,  for  an  clastic -plastic  material  subject 


to  the  yield  condition 


■  |iRT  | + wi  -  *  =  \R\ 


+  «J  S3  0 


where  and  Jg  are  the  invariants 


Ji  "  °ii 


J2  ”  2  siJ  elj 


while  the  coefficient  cr  defines  the  internal  friction  of  the  material .  The 


value  Of  is  subject  to  the  limitations  0 


<-«<-/¥ 


,  Ref.  Cl] . 


The  present  paper  will  consider  the  more  general  problem  of  a  decaying 


pressure  signal  p(x  -  Vt),  Fig.  2,  moving  also  with  superset smic  velocity  V. 


The  problem  to  be  studied  here  is,  from  a  practical  point  of  view,  more  realistic 
than  the  one  treated  in  [l],  but  the  case  shown  in  Fig.  2  is  also  considerably 


more  complex.  The  complexity  is  such  that  only  an  approximate  solution  will  be 


derived  which  is  valid  for  a  limited  distance  from  the  wove  front.  The  approxi¬ 


mation  to  be  employed  suggests  itself  when  considering  the  results  of  the  analysis 


Cl]  for  the  step  wave. 


It  was  found  in  [l]  that  the  response  changes  in  character  depending  on  the 


range  in  which  the  values  of  Poisson’s  ratio  v  and  of  the  material  parameter  ®  are 


situated,  Fig.  3.  The  values  Of  of  practical  interest  lie  between  0.1  ana  0.2, 


“I .>-«••• »  -  ,  ~cr;jsr$*rv 


[l]  Bleich,  H.H.,  Matthews,  A.T.  and  Wright,  J.P.,  Step  Load  Moving  with 

Superset srnic  Velocity  on  the  Surface  of  a  Half -Space  of  Iranular  Material , 
Tech'll'  Rpt.  AFWL-TR-65-59,  Air  Force  Systems  Command,  Kirtland  AFB, 
September  19&5* 
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•  NUMERICAL  RESULTS 

P  I  S  .^-PRINCIPAL  RANGES 


otherwise  the  value  of  the  angle  of  internal  friction  becomes  unreasonably  small, 
or  unreasonably  large.  For  values  of  Of  limited  by  0.10  <  Cf  <  0.20,  Range  III 
never  applies,  so  that  this  range  can  be  ignored,  aad  only  Ranges  I  and  II  will 
be  studied  further. 


"  SSSS£  £[ 


Consider  first  the  result  of  the  analysis  [l]  for  the  step  wave  for 
parameters  v  and  a  in  Range  II.  Defining  the  location  of  a  field  point  by  polar 
coordinates  9  and  r,  where  the  origin  of  the  coordinate  system  lies  at,  and 
moves  with  the  front  of  the  load,  it  was  found  that  the  stresses  are  functions 
of  the  angle  <p  only,  i.e.,for  a  step  load  there  is  no  dependence  on  r.  In 
Range  II, the  stress  field  shown  in  Fig.  4  was  obtained.  For  values  9  <  9, all 
stresses  vanish  while  a  stress  discontinuity  associated  with  plastic  defamations 
occurs  at  9  a  9.  The  velocity  of  this  "plasti  essure  front"  i6 


s2 


K  (1  +  20T  /?)2 

*[!♦«*  (W] 


which  defines  the  angle  9 
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(5) 


The  stress  discontinuity  at  the  front  is  such  .that  no  discontinuity  in  shear 
parallel  to  the  front  occurs.  Between  9  and  9g  there  is  no  change  in  stress 
at  all,  i .a the  mafesri n~*  **  ««■»+,  to  nlaKtic  deformations.  At  0  3  at. 

’  ~  v  *  -D 

there  is  an  elastic  discontinuous  change  in  shear,  moving  with  the  velocity  of 
elastic  shear  waves 


(6) 
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FIG.  4  CONFIGURATION  IN  RANGE  TL 


which  defines 


a  TT  - 


(7) 


In  Range  lib,  which  covers  most  of  Range  II,  see  Figs.  9b-e  of  [.l],  one  finds 
that,  from  «pg  to  a  value  ^  and  again  from  «p2  to  the  surface,  the  stresses  are 
constant  and  no  inelastic  deformations  occur.  However,  between  cp^  arid  q>? 
inelastic  deformations  do  occur,  and  the  stresses  vary  with  tp.  Details  of  the 
variation,  and  the  values  <pg  ,  and  «pg  are  given  in  Table  I  of  Ref.  [l]. 
(No  plastic  region  tp^  to  cpg  occurs  in  a  small  portion ‘of  Range  II,  called 
Range  Ila  in  [l] . ) 


As  described  above,  plastic  deformations  occur  therefore  in  at  most  two 
locations,  at  the  plastic  front  cp,  and  between  cp^  and  cp^,  .  The  prime  purpose 
of  the  preceding  description  is  to  note,  further,  that  the  plaBtic  changes  in 
stress  between  cp^  and  tp2  as  shown  in  [l]  are  small  compared  to  the  change  at 
the  front  9.  This  suggests  that  one  could  find  an  approximate  solution  that 
has  a  plastic  front  at  q>,  but  which  satisfies  everywhere  else  the  differential 
aquations  for  elastic  changes  in  stress  and  strain.  Formulating  the  problem 
in  this  fashion,  one  will  find  that  the  yield  condition  will  be  violated  in 
some  locations,  if  v  and  Of  are  in  Range  lib.  However,  one  might  accept  such  a 
solution,  if  the  degree  of  violation  of  the  yield  condition  is  not  severe. 

The  violation  can  be  recognized  by  the  excess  of  the  ratio 


n  * 


(8) 


above  unity.  Figure  5  shows  the  violation  of  the  yield  condition  for  the 
results  of  such  an  approximate  analysis  for  the  case  of  the  step  wave  fo»* 
typical  values  of  v,  Of,  and  of  the  velocity  ratio  V/cp 


It  is  seen  that  the 


violations  are  very  small  for  large  values  of  V/cp  ;  they  may  exceed  20#  for 
V/cn  <2.  As  a  first  approximation  for  engineering  purposes  one  may  accept 
violations  up  to  20#,  i.e.,  —  <  1.2.  This  order  of  violation  is  considered 

acceptable  here,  because  the  prime  purpose  of  the  present  analysis  is  to  demon¬ 
strate  the  difference  between  the  effect  of  a  step  wave,  Fig.  1,  and  the  effect 
of  a  decaying  wave,  Fig.  2.  When  results  are  obtained  later,  checks  will  be  mode 
to  confirm  the  magnitude  of  the  violation  of  the  yield  condition  in  various 
locations,  so  that  the  usefulness  and  the  range  of  applicability  of  the  approximate 
solution  can  be  judged. 


The  difference  in  the  results  for  the  step  wave  in  Ranges  I  and  II  lies  in 
the  fact  that  in  Range  I  there  is  no  plastic  shock  front,  instead  there  is  an 
elastic  shock  front  moving  with  the  velocity  of  P -waves 

4  -  !£§f  w 

The  location  of  this  front,  Fig.  6,  is  defined  by  the  angle 


«pp  *  rr  -  sin 


(1-2VM 


Following  this  shock,  the  pattern  is  in  principle  the  same  as  in  Range  II,  and 
there  is  a  division  in  Subranges  la  and  lb.  In  Subrange  lb,  which  covers  only  a 
minute  portion  of  Range  I,  see  Figs.  9b-e  of  Ref.  [l],  there  is  an  S-front  and 
a  plastic  region  of  very  small  extent  from  ^  to  ,  Fig.  6.  In  Subrange  la, 
which  covers  most  of  Range  I,  there,  is  an  S-front,  but  no  plastic  region.  In 
Subrange  la  the  solution  is  entirely  elastic  so  that  the  stress  field  without 
approximation  can  be  found  from  Ref.  [2]*^.  In  Range  lb,  the  plastic  effects  occur 


(2]  Baron,  M.L.,  Bleich,  H.H.  and  Weidlinger,  P.,  Theoretical  Studies  in 
Ground  Shock  Phenomena,  SR-19,  The  MITRE  Corporati<)n,  October  l^oT" 


FIG.  6  CONFIGURATION  FOR  RANGE  I  b 


PLASTIC 
SHOCK  FRONT 


FIG.  7  ASSUMED  CONFIGURATION 


only  in  a  minute  region  cp^  to  Tg  and  are  minor,  tf,  similar  to  the  above 
suggestion  for  Range  XI,  this  plastic  region  is  ignored,  the  solution  may  be 
obtained,  this  time  approximately,  from  Ref.  [2].  There  will  therefore  be  no 
need  to  obtain  a  new  approximate  analysis  for  Range  I  at  all. 


SECTION  II 


FORMULATION  OF  THE  APPROXIMATE  ANALYSIS 

Based  on  the  reasoning  outlined  in  Section  I,  approximate  solutions  for 
the  case  of  the  decaying  shock  wave,  Fig.  2,  will  be  formulated  for  "Range  II" 
of  the  material  parameters,  i.e.,  according  to  Eqs.  (72)  and  (111)  of  [l] 


1-2V 


-  2  *  —£=7  <  /3  «  <  1+v 

y  2(i+^)  1+v 


(n) 


This  is  achieved  by  assuming  that  plastic  deformations  occur  in  any  element  only 
at  the  time  of  the  passing  of  the  plastic  shock  front.  Fig.  7,  while  Btress 
changes  thereafter  are  elastic.  The  stresses  ahead  of  the  shock  front  vanish 
identically,  while  the  applicable  elastic  relations  between  stress  rates  and 
strain  rates  for  the  present  case  of  plane  strain  are 
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(14) 


(15) 


where  a  ,  g  ,  o  ,  T  are  the  stresses,  while  u,  v,  are  the  displacements  in 
x  y  z 

the  x  and  y  directions,  respectively.  In  addition,  the  equations  of  motion 
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must  be  satisfied.  Equations  (12,  13*  15)  may  be  integrated  with  respect  to 
time,  requiring  the  addition  of  arbitrary  functions  f^(x,y)  of  x  and  y, 
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(20) 


while  Eq.  (il)  after  integration  may  be  written 


az  "  v(<rx  +  V  *  fo*x'y* 


(SI) 


The  last  equation  is  the  ouly  relation  containing  o  ,  and  simply  defines  this 

z 

quantity.  It  is  noted  that  th?  arbitrary  functions  f^  ,  fg  and  are  not 
entirely  "arbitrary”,  but  are  related  by  a  compatibility  relation,  as  may  be 
seen  by  eliminating  u  and  v  from  Eas>  (l8),  (19)  and  (20). 

The  solution  to  the  differential  equations  (16-20)  may  be  expressed  in 
terms  of  potentials  I  and  which  satisfy  the  wave  equations 
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(22) 

(23) 


Subscripts  on  potentials  indicate  derivatives.  In  terms  of  these  potentials 
and  of  two  arbitrary  functions  g^  and  gg  of  x  and  y  the  displacements  become 

u  *  #x  “  *y +  *  .  W 


v  =  «  +  ^  +  eg(x,y) 


(25) 
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while  the  velocities  u  and  v  are 

u  »  «x  -  *y  (26) 

V  *  ^  +  I*  (27) 

the  stresses  are 

°x  *  28  fefe  '**  +  I3v  *„  -  YJ  <28> 

®y  "  20  [l-2M  *x3t  *  l-2v  *yy  +  Txy]  (2?) 

T  -  SB  ♦  *  (*«  -  V1  (30) 

with  a  given  hy  Eq.  (21). 
z 

Substitution  of  Eqs.  (24-25)  and  (28-30)  into  Eqs.  (18-20)  furnishes 
relations  between  the  three  "arbitrary"  functions  ,  fg  and  f^  ,  and  the 
two  functions  g1  and  gg  ,  i.e., 

fi’sr  >  f2-§r  >  fs  5V3r  +  «n/  <31) 

These  three  relations  and  the  Equations  (28-30)  satisfy  compatibility  identically, 
leaving  the  two  functions  g^x,y)  and  ,y)  as  arbitrary  functions  in  lieu  of 
the  related  functions  ,  fg  and  f^  . 

It  is  now  noted  that  the  velocities  u  and  v,  and  the  stresses  a%  ,  and  T 
in  terms  of  the  potentials  •  and.  i  have  exactly  the  same  form  as  in  an  elastic 
plane  problem  and  these  quantities  will  therefore  have  the  same  value  as  in  an 
entirely  elastic  problem  with  the  same  boundary  conditions.  The  plastic 
deformations  at  the  shock  front  lead,  however,  to  different  expressions  for 
u,  v  and  because  the  functions  g^  ,  gg  and  fQ  now  occur.  These  three  arbitrary 
functions  are  determinable  from  conditions  at  the  plastic  front  as  discussed  below. 
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Because  the  steady-state  problem  is  considered,  where  the  surface  pressure 


is  solely  a  function  of 
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(32) 


the  solutions  are  only  functions  of  S  and  y.  The  wave  equations  for  tne 
potentials  then  become 
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Noting  that  the  open  functions  tQ  ,  and  ,  which  are  r»wt  functions 
of  t,  can  in  the  steady-state  not  depend  on  x  either,  the  expressions  for 
displacements  and  stresses  become 


\ 


u  *  *5  '  *y  + 

(35) 

v  -  *y  +  *5  +  Sg(y) 

(36) 

*  0  [«(l  ♦  £&  SSi  -  2  *!y] 

(37) 

“y  -  o  -  2)  >%t  +  2  7{y] 

(38) 

T  -  G  [2  t;y  -  (jf  -  2)  ?;;] 

(39) 

az  -  v  (ax  *  jy)  +  ro(y) 

(40) 
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To  formulate  the  boundary  conditions,  expressions  will  be  required  for 
the  normal  and  tangential  velocities  u^  ,  u^  and  the  normal,  tangential  and 
shear  stresses  for  planes  inclined  at  ai<y  angle  tp  (see  Fig.  8) 

6 

V"  *  8in  *  “  C0»  <?  9|y  "  C0B  *  ?{{  -  ¥  ?gy  ( kl ) 

A 

7“  *  -  cos  p  #gg  -  sin  p  »|y  -  sin  p  tgg  +  cos  p  7gy  (1»2) 

(1  -  ~~  sin2  p)  -  2  coe  2p]  -  2  sin  2*  #?y  + 

+  2  cos  2*  ¥gy  +  (1^  -  2)  sin  2?  ?gg  (43) 

ff 

T  -  P  i.om  p 

(1  -  cos  f )  +  2  cos  2p]  « gg  +  J  Bin  Sp  - 

-  a  cm  ap  »{y  -  (i|  -  2)  sin  2p  hi  W 

TN  2 

-  (IT  -  2)  sin  2?  »gg  +  2  cos  2*  #gy  +  2  sin  Sp  fgy  - 

-  <i(j  -  a)  cm  a?  T{{  (45) 

To  state  the  boundcry  conditions  of  the  problem  properly,  and  to  be  able 
to  make  statements  on  the  character  of  the  solutions,  the  differential  equations 
for  i  and  t  must  be  discussed.  While  both  differential  equations  are  hyperbolic, 
the  characteristics  for  *  are  steeper  than  the  boundary  p  defining  the  plastic 
front,  <  p  .  The  solution  for  f  in  the  region  of  interest,  (p  <  p  <  it,  Fig.  7, 
is  therefore  defined  by  statements  on  the  boundaries  cp  «  9  and  9  »  n  .  §  ,ill  be 
a  continuous  function  if  the  prescribed  condition?  at  the  boundaries  are  con¬ 
tinuous.  This  being  the  case,  it  will  be  possible  to  use  an  expansion  in  a 
Taylor  series  in  {  and  y  for  a  region  near  the  origin  5  *  y  ®  0. 
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On  the  other  hand,  the  characteristic  directions  for  the  differential 
equation  on  f  are  flatter  than  the  boundary,  90  >  9  .  The  differential  equation 
does  therefore  not  apply  at  9  »  9g  ,  and  the  solutions  Y  for  9  >  9g  »  and  Y  for 
9  <  9  <  are  not  fully  defined  by  prescriptions  on  9  ■  t?  and  9  «  9  .  The  two 
functions  Y  and  Y  cast,  however,  lead  to  values  of  stress  and  velocities  that 
satisfy  momentum  considerations,  which  furnishes  additional  conditions  for  their 
determination.  While  at  least  some  of  the  derivatives  of  Y  and  Y  at  9  «  9  will 
not  be  continuous,  each  of  these  potentials  in  its  region  will  be  continuous, 
and  each  of  the  potentials  can  again  be  expanded  in  a  Taylor  series. 

The  functions  $,  Y  and  Y  are  then  subject  to  the  following  conditions.  On 
the  surface  the  normal  and  tangential  stresses  are  prescribed.  Thus  for  an 
exponentially  decoying  surface  load, 


at  y  »  0,  j  <  0 

fry  »  -  PQ  ew5  (it  >  0)  (4S) 

T  *  0  (47) 

At  the  plastic  front  9  »  9  conservation  of  momentum  requires 


frR  »  p  c  Ujj  (48) 

«  0  (49) 

where  c,  9  are  given  by  Eqs.  (4)  and  (5),  respectively.  In  addition  the  yield 
condition  most  be  satisfied  at  the  plastic  front  so  that  the  stresses  and 
are  related  to 
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at  tp  «  q> 

■  "  (rils^s)  au  (50> 

*  (H»^)  ®k  (51) 

Due  to  the  fact  that  different  egressions  ¥  and  f  for  the  solution  of 
the  second  of  Eqs.  (33)  vill  be  used,  it  is  necessary  to  consider  additional 
conditions  at  the  location  «p  «*  .  Momentum  considerations  require  that 

discontinuities  in  the  shear  stress  and  the  tangential  component  of  the 
velocity  u^  are  proportional.  Using  the  symbols  ,  etc.  for  the  discontinuities 


at  <p  »  <pg 


4tk  -  e  c3  4^ 


while  all  the  normal  and  tangential  stresses  and  the  normal  velocity  at 
«p  =  cpg  must  be  continuous 

at  (p  »  <ps 

A<?n  «  0  (53) 

a  0  (5*0 

«  0  (55) 

AUj^  »•  0  (5^) 

The  above  equations  complete  the  available  conditions  for  the  determination 

_  as 

of  the  potentials  ¥  anil  ¥.  Aftor  th*?  potentials  ere  obtained  the  arbitrary 
function  fQ(y)  appearing  in  3q.  (U&)  can  be  obtained  from  Eq.  (51)*  The 
arbitrary  functions  g^(y)  and  gg(y)  in  E<5S  •  (35 )  and  (36)  are  defined  by  the 
fact  that  the  displacements  must  vanish  at  the  plastic  front 


u  a  7  *  0 


IS 


M 


'  W?: 


SECTION  III 


SOLUTION  OF  THE  DIFFERENTIAL  EQUATIONS  BY  TAYLCR  SERIES. 


The  differential  equations  for  the  potentials  to  be  solved  are 


»„ '  ^  -  1)  {t£  '  0 

V  -  - »'« - 0 


The  boundary  conditions,  Eqs.  (46-47)  on  the  surface  give 


»  o,  5  <  o 


<4  -  2> ’«  +  2  ?!y  -  -  T 

'  i  ! 


2  »{y  -  <l|-  2)»S{-  0 


While  the  conditions,  Eqs.  (48-50)  at  the  plastic  front  give  three  equations 


sin  2(p  #5y  + 


[>!  -  2  +  .in2  f  (4  -  2M2  -  §K  +  (~g~)  1 
+  [2  +  (f^  -  4)  sin2  $)  ?5y  +  (^r-“)  sin  2?  ?, 


(tf*  -  2)  sin  2?  *g5  +  2  cos  +  2  sin  2tp  Tjy  -  (1^  ”2)  cos  2(p  *  0  (63) 

S[co»  2$  +  <l|  -  t?)  (h  *  f$)  -  |  +  £  .In2  f]  #„  + 


+  2  sin  29  -  2  cos  •V-'sfi  -2)  sin  2a  fgm  u  0 


‘•<*  w  '  -t 


Finally,  the  five  conditions,  Eqs.  (52-56)  at  the  S-front  lead  to  only  one 
condition 


at  y  -  5  taa  «pg 

!-»«)*  *“  *s  (?{y-V-o 

The  surface  load  is  now  expanded  in  a  series  in  £ 

VU'  -  V  £ 

m»0 


(65) 


(66) 


applicable  for  5  <  0,  while  the  potentials  #,  Y  and  Y  are  expanded  in  double 
series 


The  differential  equations  and  the  conditions,  Eqs,  (58-65)  are  homogeneous, 

_  as 

i.e.,  they  contain  consistently  only  second  derivatives  of  #,  Y  and  Y  with 
respect  to  ?  and  y,  so  that  the  recurrence  equations  for  the  values  a,  b,  c 
couple  coefficients  only  for  which  i+j  has  the  same  value.  In  view  of  Eq.  (66), 
these  terms  couple  only  with  the  term  m  »  i+J-2  of  the  expression  for  the 
surface  load.  The  potentials  may  therefore  be  written  as  a  series  of  polynomials 


Z* 

m«0 


,(«) 


i(») 


m+2 


'-2> 

m«0 


w(n) 


i-0 

m+2 


I'?'  «* 


.  m+2-i 


,  m+2-i 


■I* 

m»0 


*(m) 


i-0 

m+2 


«  £c(»)  5s-  yB+2-i 

i»0 


(68) 


(69) 


(70) 


«•  ■»« 

« 

Y 

=  7 

iW 

(67) 

. 

f 

? 

_Ci ,S_ 

^  | 

J 
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The  case  m  *  0  corresponds  to  the  step  load  and,  while  quite  simple,  is 
treated  differently,  see  Appendix  I.  For  values  a  >  1  there  are  (3m  +  9) 
unknown  coefficients.  Substituting  the  second  derivatives  of  Eqs.  (68-70) 
into  Eqs.  (60-65}  yields  six  equations: 


<4  ■  s)i”  +  2)(a  +  1J  a«*2  +  2(0  +  X)  Vl  •  -  l  »! 

(7D 

2(»  *  1)  -  (>^  -  2)(m  +  2)(m  +  1)  o^g  -  0 

(72) 

a+2 

X(bJ  ”  CJ}  J  (tan  fs)^2"*3  “  0 

(73) 

J-l 


(m  +  2)(a  +  1)  +  N^bJ  +  (m  +  l)(tan  <p)m  [n^  +  N^J 

rr  '  1 

+  L  (m  +  3  “  J,(t*n  ’),J  1  {  +  2  -  N1  +  (*}  -  i)  cot  *  N2J  a 


m+3-J 


+  [(m  +  2  -  j)  +  (J  -  1)  cot  9  hJ  V3-}}  "  0  <7*0 

(a  +  2)<a  +  1)  [N5am+2  +  Ngb^J  +  (a  +  l)(tan  if  [l^  +  + 

m+l 

+  ^(i  +  3  -  J)(tan  i)*'1 1  [(m  +  2  -  J)  N5  +  (J--  1)  cot  $  N6J  a^^  + 


(a  +  2  -  J )  Nq  +  (j  -  1}  cot  $  ^  bm+3-j| 


(75) 
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(B  +  2)(b  + 1)  [yH+2  +  s10*Bt2J  *  (»  +  «(«“  if  [vi  -  %bi]  + 

m+i  r  r  i 

+  £{a  +  3  -  J)(taa  $)*ml  {  [(m  +  2  -  j)  N9  +  (j  -  1)  ^  cot  $j  + 

j-2 

+  [(n  +  2  -  J)  N10  -  0  -  1)  N6  cot  9]  bm+3_J  «  0  (76) 


where 

W1  *  Kg  -  2  +  8in2  9  (*  * 

“2  *-(rr-)Blna* 

N3  -  2  +  (J^  -  4)  sin2  9 

“it  ’  (~~2  ~)  "ln  ^ 

1?5  -  (Np  -  2)  sin  29 

Ng  «  2  cos  29 

Kj  m  2  sin  29 

“8  -  -  0|  -  2)  COS  29 

N9  »  2[cos  29  +  -  &)  (HIt0  -  I  +  &  ®in2  9] 

N10  *  '  "  2)  8i0  2$ 


Substituting  Eqs.  (68-70)  into  the  wave  equations,  and  equating  coefficients 
of  like  powers  of  S  and  y,  yields  (3m  +  3)  recurrence  equations 
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aa-2q+2  * 


rf-^wTT  Pm,l 

4*1 


am-2q+l  "  ^  ”  1)<3  Vl  *2q-2,a  TT  ^ 


Ws  "  <*€  “ 1)9 


bm+2  TT  Pn>4 


Wl  “  *  1)<3  Vl  4 


TT  § 

2q-2,a  j  j  m,4 


m-2q+2 


-  i'q  w  IT  pm 


cm-2q+l  "  ^  '  1*<I  cm+l  ®aj-2,m  TT  Pm,* 

4*1 


where 


a, 4  *  '"3(54 

^T) 

1 

L 

a  (a  +  2  -  24) 

(*  + 

3  -  3*, 

l 

1m,4  24(24 

+  1) 

i 

1 

if 

2q  - 

2  f  m 

®2q-2,a  “  0 

if 

2q  - 

2  *  m 

1  <  q  < 

mm  *•  mm 


2a  +  3  +  (-  lj 
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Equations  (71-76)  and  Eqs.  (78)  give  the  required  number  of  (3a  +  9) 
equations  for  the  (3a  +  9)  unknown  coefficients  a^  ,  b^  and  . 

Substitution  of  Eqs.  (78)  into  Eqs.  (?1)  to  (76)  reduces  the  system  to 
a  set  of  six  linear,  nonhomogene ous,  simultaneous  equations,  involving  only 
the  unknowns  9^  ,  ,  ba+1  ,  c^  ,  .  After  solving  for 

these  values,  the  stresses,  velocities  and  displacements  fov  each  of  the 
m  components  of  the  surface  load  may  be  computed  at  any  location  (?,y) 
behind  the  front  at  9 .  As  an  illustration, the  details  for  the  case  m  *  1 
are  written  out  in  full  in  Appendix  37. 
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SECTION  IV 


HUSflSRIGAL  RESULTS  AND  DISCUSSION. 

Numerical  results  fro®  the  relations  derived  in  the  previous  section  were 
obtained  on  a  digital  coc^uter*  truncating  the  series  solution  by  retaining  only 
the  terns  up  to,  and  including  m  ■  5.  The  stresses  found  are  presented  in 
Table  IX  -  XII  and  in  the  accos^enying  figures,  for  the  values  V/'c .  »  2  and  $  for 
conciliations  of  the  par  alters  ©  ■  0.10,  0.20,  V  «  0*12$,  0.2$,  0.333*  In 
addition  the  case  v/cp  »  2,  ar  *  v  ■  0.0$,  which  lies  very  close  to  the  linlt 
between  Ranges  II  and  III,  is  also  given*  The  caabinations  of  parameters  for 
which  results  are  available  are  marked  in  Fig*  3* 

The  neglect  of  a  plastic  region,  which  may  occur  between  the  S -front  and 
the  surface,  restricts  the  applicability  of  the  resulting  approximate  solutions 
to  distances  from  the  front  of  the  applied  load  where  the  ratio  n,  Eq.  (8), 
exceeds  unity  only  moderately.  To  check  this  matter  the  Tables  Indicate  the 
value  n  for  all  points  where  stresses  arc  listed.  However,  the  truncation  of 
the  series  further  limits  the  validity  of  the  results  to  ||sC|  <  l.$.  The 
truncation  at  value  m  *  $  became  necessary,  because  for  larger  values  accuracy 
troubles  developed  in  the  computations.  The  numerical  results  presented  contain, 
for  each  combination  of  parameters,  only  the  range  where  the  approximate  solution 
is  meaningful,  allowing  for  both  approximations. 

1.  Discussion  of  typical  results  for  V/ep  m  2 

The  case  a  ■  0.10,  v  °  0.12$,  V/cp  «  2.0  is  typical  and  is  therefore  selected 
for  discussion.  Figure  9  shows  the  nondimensional  history  of  the  principal  stress 
cr^  on  the  surface  y  *  0,  and  for  two  values  of  the  depth,  (Ay  »  0.0$  end  0.40.  The 
values  ere  plotted  against  tan  f  *  (pit  taa  v] ,  where  t  is  the  time  elapsed 
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RESUL 


after  arrival  of  the  surface  load  above  the  target  point*  The  curve  for  y  *  0 
shows  Just  the  history  of  the  applied  load,  while  the  other  curves  indicate  a 
discontinuous  increase  in  (-  0^)  at  the  arrival  of  the  plastic  front,  followed 
by  a  decrease  in  (-  0^).  This  *s  subsequently  followed  by  a  second  discontinuous 
increase  in  (-  er^)  at  arrival  of  the  S-front,  followed  by  a  gradual  decrease. 
Attention  is  drawn  to  an  import  ant  difference  between  the  ciurves  for  |iy  »  0.05 
and  |*y  »  O.hO.  The  fonaer  is  representative  of  the  situation  at  small  depths, 
where  the  maximum  value  of  (-  a^)  occurs  at  the  second  peak  due  to  the  passing 
of  the  S-front.  The  curve  for  |iy  *  O.hO  is  representative  of  the  situation  at 
larger  depths,  when  the  first  peak  due  to  the  passing  of  the  plastic  front  gives 
the  maximum  value  of  (-  0^).  The  equivalent  occurs  in  the  case  of  a  decaying 
load  on  an  entirely  elastic  medium,  and  Jig.  9  shows  the  comparative  histories 
of  (-  e^)  at  the  two  depths  |iy  <*  0.05  and  0.40  in  dashed  lines.  The  arrival 
time  of  the  first  peak  in  the  elastic  cases  is  of  course  earlier,  but  the  relative 
values  of  the  first  and  second  peaks  bear  the  same  relation  as  for  the  elastic- 
plastic  material. 

Further  details  for  the  above  values  of  tbe  parameters  are  given  in  Table  II 
and  its  acconpanyiag  figure.  Table  H  shows  the  values  of  the  principal  stresses 
°1  ’  °2  *n  ^  P^-ha®  38  v®!1  38  third  principal  stress  ,  and  the 
angle  8  between  the  direction  of  and  the  {  axis.  In  addition,  the  values 
of  n,  indicating  the  violation  of  the  yield  condition,  are  listed.  The  in¬ 
formation  In  Table  S  is  given  for  a  rectangular  network  of  points,  1-10,  selected 
such  that  one  point  at  each  of  the  depth  levels  falls  on  the  plastic  front. 
Additional  points,  11  to  13,  lie  on  the  S -float.  At  point  I  two  values  for  eacn 
quantity  are  given,  because  the  stress  field  has  a  singularity,  i.e.  the  values 
depend  on  the  direction  of  approach  <p,  and  the  results  differ  depending  on 
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whether  9  ^  9g  .  (The  two  sets  of  values  at  point  1  immediately  behind  the 
frjnt  of  the  applied  load  are  of  course  equal  to  those  due  to  a  step  load  with¬ 
out  decay. ) 

The  values  n  in  the  last  coluir  of  Table  2  permit  evaluation  of  the 
approximation  achieved.  There  is  a  moderate  violation  of  the  yield  condition 
at  point  1,  n  *  1.19*  As  mentioned  above,  this  point  represents  the  values  due 
to  a  step  load  without  decay,  for  which  case  results  are  known  without  approximation. 
To  judge  the  seriousness  of  the  violation  expressed  by  n  =  1J.9  >  1  Table  1  gives 
the  exact  values  of  the  stresses  and  angles  8.  It  is  seen  that  the  differences 


TABLE  I  Exact  values  at  {  *  y  »  0,  for 

V  ■  0.125,  «  «  0.10,  V/cp  *  2.0 


Location 

-  cr,/p 
X'^o 

-°A 

jaP 

1 

e 

9  <  9  <  9g 

0.94U 

O.58O 

0.580 

65.03° 

9  a  ff 

1.00 

0.598 

0.632 

90.00° 

in  the  principal  stresses  are  not  large,  the  maximum  being  9$*  The  errors  at  the 
shock  front  are  very  small,  less  than  2 j>.  While  the  values  8  listed  in  Table  1 
agree  with  those  in  Table  2,  it  is  noted  that  in  the  exact  case  the  value  8  varies 
between  (p  »  cpg  and  9  a  IT  from  76.76°  to  90°,  while  the  approximate  value  for  9 
in  this  entire  range  is  90°.  While  this  difference  in  direction  may  not  be  of 
great  practical  importance,  it  is  percentage -wise  of  the  order  of  tne  error  in  n, 
and  about  twice  as  large  as  the  error  in  the  stresses. 

Moving  away  from  point  1  at  the  origin  to  adjacent  points  2,  5,  H  at 
H§  =*  -  0.2X5,  and  to  further  points  at  wS  =  -  0.859  and  =  -  I.289  it  is  noted 
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TABLE  III 


Results  for  v  =  0.125,  ©  *  0.10,  V/cp  =  5.0 


IS  4  3  2  I 


Point 

■  ' 

-  5*5 

ity 

-  °iAo 

-  *A> 

-  a3/po 

9° 

a 

1  * 

0 

0 

0.9938 

0.6103 

0,6103 

80.20 

1.00 

1  ## 

0 

0 

1.00 

0.6041 

0.6103 

1.02 

2 

0.292 

0 

0.7469 

0.5621 

0.5734 

90.00 

0.550 

3 

O.876 

0 

0.4161 

0.5018 

0.5245 

90.00 

0.396 

4 

1.167 

0 

0.3082 

0.4794 

0.5082 

90.00 

0.835 

5 

0.292 

0.050 

0.9523 

0.5848 

0.5848 

80.28 

1.00 

6 

0,676 

0.050 

0.5456 

0.4974 

0.5231 

55.31 

0.155 

7 

1.167 

0.050 

0.4942 

0.3878 

O.5029 

104.1 

0.463 

8 

O.876 

0.150 

0.8751 

0.5374 

0.5374 

80.28 

1.00 

9 

1.167 

0.150 

O.6636 

0.4924 

0.5053 

71.80 

0.574 

10 

1.167 

0.200 

0.8391 

0.5152 

0.5152 

80.28 

1.00 

11*** 

0.292 

0.039 

0.9016 

0.5771 

0.5813 

87.86 

0.904 

12*** 

0.876 

:  0.116 

0.7393 

0.5240 

0.5294 

79.77 

0.685 

13*** 

1.167 

0.154 

0.6756 

0.4944 

O.5058 

72.94 

0.606 

14*** 

1.459 

0.193 

0.6250 

0.4602 

0.4837 

64.46 

0.568 

15 

1.459 

Q 

0.2213 

O.4599 

O.4949 

yo.oo 

1.27 

TABLE  IV 


Results  for  v  »  0.125,  ®  *  0.2Q,,  V/cp  =»  2.0 


Point 

-  v& 

uy 

-  Vpo 

-«A 

1° 

n 

1  *' 

0 

0 

0.9403 

0.3631 

0.3631 

61.63 

1.00 

1  *# 

0 

0 

1.000 

0.3031 

O.363I 

90.00 

l.afi 

2 

0.278 

0 

0,7575 

0.2869 

0.3307 

90.00 

0.945 

3 

0.556 

0 

0.5737 

0.2739 

U.;t06l 

90.00 

0.713 

4 

0,926 

0 

0.3954 

0.2607 

0.12821 

90.00 

0.386 

5 

0.278 

0.15 

0.9241 

0.3568 

0.3568 

61.63 

1.00 

6 

0.556 

0.15 

0.7841 

0.2424 

0, 3250 

94.82 

1.08 

7 

0.926 

0.15 

0.5412 

0.2386 

0.2942 

95.80 

0.750 

8 

0.556 

0.30 

0.9083 

0.3507 

0.3507 

61.63 

1.00 

9 

0.926 

0.30 

0.7516 

0.1883 

0.3108 

99.19 

1.18 

10 

0.926 

0.50 

0.8671 

0.3348 

0.3348 

61.63 

1.00 

11*#* 

0.278 

0.096 

0.9239 

0.2650 

0.3465 

92.98 

1.17 

12##* 

0.556 

0.192 

0.8590 

0.2272 

0.3315 

95.90 

1.19 

13*^ 

0.926 

0.321 

0.7873 

0.1782 

0,3135 

99. 58 

1.25 

ill 

LI _ 1 

1.30 

0 

0.2669 

0.2505 
L_ _ J 

0.2648 

90.00 

0.057 

*  Values  at  9  <  <  <ps  **  Values  at  ipg  <  cp  <  n 

*#*  Values  at  9  ■  9*  .  At  9  ■  9g  ,  At  *  ♦  0.403  p  . 


Points 

-u5 

)ij 

I 

Q 

0* 

- 

'  “3/p° 

n 

0 

0 

0.9915 

O.3828 

0.3828 

79.0i» 

■  1 

0 

0 

1.000 

0.3743 

0.3828 

90,00 

2 

0.258 

0 

0.7724 

0.3432 

0.3505 

90.00 

Ifgi!  I 

3 

0.775 

0 

0.4605 

O.2999 

0.3061 

90.00 

■ 

4 

1.30 

0 

0.2669 

O.C720 

0.2784 

90.00 

5 

0.258 

0.05 

0.9786 

0.3773 

0.3778 

79*04 

1. 00 

6 

0.775 

0.05 

0,5912 

0.3125 

0.3212 

91.61 

0.647 

7 

1.30 

0.05 

0.3475 

0.2799 

0.2867 

91*97 

0,204 

8 

0.775 

0.15 

C.9533 

O.368I 

0.3681 

79*04 

i.00 

9 

1.30 

0.15 

0,5794 

0.2956 

0.3123 

95-74 

0.671 

10 

1.30 

0.252 

0.9283 

0.35^ 

0.3584 

79.0* 

1.00 

11*** 

0.258 

0.034 

0.9159 

0.3574 

0,3683 

91.05 

0.973 

12*** 

0.775 

0.10 

0.7611 

0,3258 

0.3414 

93*4? 

0.864 

13*** 

1.30 

0.6500 

0,2985 

96.58 

0.776 

14*** 

1.55 

0.5983 

0,2855 

0.3105 

98.39 

0.728 

15 

1»55 

0.1966 

0.2616 

0.2683 

90.00 

0.272 

*  Values  at  9  <  «p  <  5(»g  .  **  Values  at  !pg  <  <p  <  ff 

***  Values  at  9  =  9q  .  At  9  «  9e  ,  At  *  +  0.178  . 
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Point 

- 

»*y 

-  V»o 

ma 

5° 

n 

1  * 

0 

0 

0.9791 

0.6012 

0.6012 

62.10 

1.00 

1  ** 

0 

0 

0.5803 

0.6012 

90.00 

1.08 

2 

0.283 

0 

0.7533 

c.5033 

0.5203 

90.00 

0.786 

•V 

m, 

U 

0.756 

0 

0.4695 

0.4135 

0.4269 

90.00 

0.224 

** 

1.133 

0 

0.3195 

0.3648 

0.3772 

90.00 

0.286 

5 

0.283 

0.15 

0.9566 

0.5874 

0.5874 

62.10 

1.00 

6 

0.756 

0.15 

0.6444 

0.4170 

0.4666 

100.6 

0.782 

7 

1.133 

0.15 

9.4450 

0.3629 

0.4034 

110.4 

0.339 

8 

0.756 

0.40 

0.9203 

0.5651 

0.5651 

62.10 

1.00 

9 

1.133 

o.4o 

0.6332 

0.4645 

0.4682 

63.69 

0.615 

10 

1.133 

0.60 

0.8923 

0.5479 

0.5479 

62.10 

1.00 

11*** 

0.283 

0.085 

0.8973 

0.5126 

0.5559 

94.91 

1.C7 

12*** 

0.756 

0.228 

0.7668 

0.4072 

0.4925 

102.8 

1.13 

13*** 

1.133 

0.342 

0.6885 

0.3299 

0.4501 

108.3 

1.24 

14 

1.30 

. . «— . 

0  . 

L, _ 

0.2669 

_ 

0.3474 

0.3597 

90.00 

0.516 

*  Values  at  9  <  (?  <  <pg  **  Values  at  <pg  <  <p  <  n 

***  Values  at  9  «  tpg  .  At  <p  =  qpg  ,  At  =  ♦  0.437  pQ  . 
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TABLE  VII 


Results  for  v  =»  0.25,  »  =  0.10,  V/cp  =  f  >0 


15  4  3  2  1 


Point 

-  t*S 

i*y 

-  eA 

-  V»o 

-  °3/po 

gO 

n 

0 

0 

0.9968 

0.6121 

0.6121 

79*23 

1.00 

0 

0 

1.000 

0.6089 

0.6121 

90.00 

1.01 

2 

0.262 

0 

0.7692 

0.5313 

0.5350 

90.00 

0.742 

3 

0.787 

0 

0.4543 

0,4245 

0.4297 

90.00 

0.124 

4 

1.312 

0 

0.2633 

0.3580 

0.3652 

90.00 

0.576 

5 

0.262 

0.05 

0.9795 

0.6015 

0.6015 

79.21 

1.00 

6 

0.787 

0.05 

0.5829 

0.4514 

0.46 73 

92.70 

0.453 

7 

1.312 

0  ^5 

0.3810 

0.3421 

0.3870 

90.42 

0.220 

8 

0.787 

0.15 

0.9437 

0.5807 

0.5807 

79.21 

1.00 

9 

1.312 

0.15 

0.5723 

0.4360 

0.4512 

101.0 

0.512 

10 

1.312 

0.25 

0.9131 

0.5607 

0.5607 

79.21 

1.00 

11*#* 

0.262 

0.031 

0.8943 

O.56Y8 

0.5732 

91.39 

0.919 

*1  OJLMJ; 

0.787 

0.092 

0.7181 

0.4970 

0.5070 

95.12 

0.725 

1.312 

0.153 

0.5794 

0.4374 

0.4531 

101.2 

0.530 

1*575 

0.183 

0.5226 

0.4106 

0.4301 

105,6 

0.439 

15 

1.575 

0 

0.1899 

0.3321 

0.3404 

90.00 

0.961 

*  Values  at  f  <  <p  <  9g  .  **  Values  at  9g  <  <p  <  n 

***  Values  at  9  =  9g  .  At  9  «pg  ,  A?  =  +  O.I85  pQ  » 
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TABLE  VIII 


FIG.  16  PRINCIPAL  STRESS  6\  FOR 
V  *0.25,  CX.-0.20,  V/ip*2.0 


42 


»  168.47 


Point  - 


IX*** 

12**» 

13*## 


0.246 

0.737 

1.23 

0.246 

0.737 

1.23 

0.737 

1.23 

1.23 

0.246 

0.737 

1.230 

1.4/4 

1.474 


0.05 

0.05 

O.05 

0.15 

0.15 

0.25 

0.029 

0.086 

0.143 

0.171 


0.9942 

O.6085 

0.9930 

O.9021 

0.7384 

0.608b 

0.5553 

0.2174 


0.3846 

0.3808 

0.3094 

0.2097 

0.1476 

0.3844 

0.2488 

0.1723 

0.3839 

0.2552 

0.3834 

0.3442 

0.2803 

0.2258 

0.2014 

0.1241 


-  Vpo 


0.3846 

0.3846 

0.3123 

0.2115 

0.1485 

0.3844 

0.2556 

0.1765 

0.3839 

0.2553 

0.3834 

O.3510 

0,2941 

0.2480 

0.2285 

0.1248 


78.49 
9?.  52 
92.8o 
78.49 
78.49 
"8.49 
91.34 
94.42 
y8.nf 
loo.] 
90.00 


O.969 

O.896 

0.695 

1.00 

0.937 

A.  00 

0.911 


0.993 
0.993 
1 .00 
O.V75 


Values  at  cp  <  9  <  #*  vajUf?K  at  ep„  <  <  n 

S>  •"* 

***  Values  at  <p  *  ,  At  Cp  »  ^  *  +  O.I58  pQ  . 
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AT  SHOCK  FRONT 


L_4_l _ 1 _ L _ 1 _ ( _ 1 _ I _ L 

0.5  1.0 


FIG.  17  PRINCIPAL  STRESS  <g,FOR 
V  *0.25.  (X  *  0.20;  y/  *5.© 


<p  =>  150.42° 


Point 

~ 

-  »s 

uy 

-  °iK 

- 

-  °3/po 

e° 

n 

1  * 

0 

0 

1.002 

0.6155 

0.6155 

60.42 

1.00 

1  ## 

0 

0 

1.000 

0.6179 

0.6155 

90.00 

0.991 

2 

0.264 

0 

O.7678 

0.4998 

O.4987 

90.00 

0.878 

3 

0.705 

0 

0.4941 

0.3605 

O.36U 

90.00 

0.633 

4 

1.057 

0 

0.3458 

0.2849 

O.2865 

90.00 

0.378 

5 

0.264. 

0.15 

0.9978 

0.6127 

0.6127 

60.42 

6 

0.705 

0.15 

0.6797 

0.3957 

0.4343 

100.8 

1.02 

7 

1.057 

0.15 

0.4788 

0.3089 

0.3385 

102.8 

O.io6 

8 

0.705 

o.4o 

0.9902 

0.6081 

O.6081 

60.42 

i.OG 

9 

1.057 

0.40 

0.6953 

0.4602 

0.4605 

60.65 

0.339 

10 

1.057 

0.60 

0.9842 

0.6044 

0.6044 

60. 42 

1.00 

11**# 

0.264 

0.068 

0.8827 

0.5298 

O.5469 

94.79 

1.02 

12*** 

0.705 

0.182 

0.7310 

0.3992 

0.4526 

102.4 

1.13 

13*** 

1.057 

0.273 

0.6403 

0.3066 

0.3919 

107. e. 

1.29 

14 

1.30 

0 

0.26&9 

0.2447 

0.2468 

90.00 

J.161 

_ 

u  ■ 

,  .  1 

.  _ 

..  .  ....  — j 

U - ,-J 

_ 

*  Values  at  ^  <  (j)  <  cpg  **  Values  at  «pg  <  <p  <  TT 

***  Values  at  (p  =  .  At  9  »  «pg  ,  AT  =  4-  0.392  pQ  . 
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TABLE  XI 

Results  for  v  *  0.333,  *  =  0.10,  v/e 


Point 


0 
0 

0.2 

0.7 

1.2 
0.2 
0.7j 

1.2 

0.7^5 

1.2 

1.2 
0.2l 
0.7 1 
1.2* 
l.W 


0 

0 

0 

0 

0 

0.05 

0.05 

0.05 

0.15 

0.15 

0.25 

0.025 

0.075 

0.125 

0.15 


-  Vpo 

-  V»o 

•  VPo 

f  9° 

n 

0.9999 

0.6l4o 

0.6l4o 

78.61 

1.00 

1.000 

0.6139 

0.6l4o 

90.00 

1.00 

0.7802 

0.5038 

0.5040 

90.00 

0.892 

0.4746 

0.3507 

0.3512 

90.00 

0.607 

0.2847 

0.2555 

0.2561 

90,00 

0.210 

0.9965 

0.6119 

0.6119 

78.61 

1.00 

0.6112 

0.4118 

0.4168 

94.68 

0.790 

0.3711 

0.2937 

0.2074 

97*20 

0.454 

0.9897 

0.6078 

0.6078 

76*61 

1.00 

0,6025 

0.4134 

0.4139 

78.31 

0.762 

0.9830 

0.6037 

0.6037 

78.61 

1.00 

0.8641 

0.5517 

0.5545 

91*93 

0.960 

0.6946 

0.4460 

0.4559 

96.88 

0.882 

0.5515 

0.3593 

0.3790 

103.5 

0.820 

0.4945 

0.3212 

0.3472 

107.4 

j.dou 

u.tfiy?  j 

O.2130 

0.2202 

90. 00 

0,061 

Values  at  c?  <  9  <  cpg 


**  Values  at  cps  <  tp  <  n 


***  Values  at  <p  *  fp*  At  9  >  f „  ,  Av  =  +  0.161  p  . 


•AT  SHOCK  FRONT  $  «  $ 


FIG.  19  PRINCIPAL  STRESS  6 ,  FOR 

*9*0.3333,  C<*O.IO,  tycp*5.0 


TABLE  XT! 

Stesults  for  v  =  0.05,  0 t  »  Q.05,  v/cp  «  2.0 


Point 

-1*5 

lAy 

-  V»o 

-  °A> 

-  °3/po 

e° 

n 

1  * 

0 

0 

0.9668 

■HH 

0.75?7 

69.33 

1.00 

1  #* 

0 

0 

1.000 

0,7527 

90.00 

1.24 

2 

0.133 

0 

0.8759 

0.7461 

90.00 

0.741 

3 

0.398 

0 

0.6953 

0.6719 

0.7351 

90.00 

0.304 

4 

0.663 

0 

0.6775 

0.5154 

0.7264 

90.0c 

1.151 

5 

0.133 

0.05 

0.9338 

0.7270 

0.7270 

69.33 

1.00 

6 

0.398 

0.05 

0.7420 

0.6810 

0.7151 

88.79 

0.286 

7 

0.663 

0.05 

0.6668 

0.5687 

0.7058 

92.22 

0.728 

8 

0.398 

0.15 

0.8721 

0.6790 

0.6790 

69.33 

1.00 

9 

0.663 

0.15 

O.6952 

0.6381 

0.6681 

86.44 

0.286 

10 

0.663 

0.25 

O.8156 

0.6350 

0.6350 

69.33 

1.00 

11*** 

0.133 

0.049 

0.9640 

0.6937 

0.7275 

90.43 

1.24 

12*** 

0.398 

0.146 

O.8989 

0.6453 

0.6804 

91.30 

1.24 

13*** 

0.663 

0,243 

O.8396 

0.6011 

0.6372 

92.15 

1.24 

14 

0.795 

0 

0.6683 

0.4512 

0.7227 

90.00 

1.56 

*  Values  at  $  <  ?  <  <|>s  **  values  at  <ps  <  c?  <  n 

***  Valr-s  at  (p  =  .  At  «p  =  tpg  >  At  ■  ♦  0,556  pQ  . 
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that  the  largest  violations  n  >  1  occur  at  points  on  the  S-front.  Por  the 
values  of  l»C  considered  above  the  violations  are  not  greater  than  at  point  1. 

The  results  for  these  three  values  of  may  therefore  be  accepted  as  being  in 
a  region  of  good  approximation.  This  desirable  situation  does,  however,  not 
continue  for  larger  values  of  (-  )•  Points  13  and  14,  both,  show  rapid  increases 

in  n  indicating  poor  approximations  for  j  j  >  1.30.  In  the  case  dlscuseed,  the 
violations  on  the  S-front  are  more  severe  than  on  the  surface.  This  is  typical 
for  V/cp  ■  2.  For  much  larger  values  of  this  parameter,  e.g.,  V/cp  «  5,  Table  III, 
it  is  seen  that  the  violations  on  the  surface  y  »  0  control  the  range  of  validity 
of  the  results. 

Having  established  that  the  solution  in  the  range  of  the  variables  and  y 
covered  in  Table  II  is  a  reasonably  good  approximation,  the  (approximate)  solution 
found  for  the  decaying  pressure  pulse  can  be  compared  with  the  (exact)  solution  [l] 
for  the  step  wave  to  assess  the  importance  of  the  decay.  For  the  step  pulse  the 
peak  stress  (-  0^)  at  each  target  point  occurs  at  the  arrival  of  the  S-front,  the 
peak  value  remaining  at  0^/po  »  -  1  for  all  depths.  Figure  10  shows  the  decrease 
of  the  principal  stress  with  depth  for  the  decaying  pulse.  It  is  seen  that 
(-  <?,/p_)  in  the  range  shown  decreases  only  to  about  75$*  It  is  important  that 
this  very  moderate  decrease  is  not  Just  due  to  the  fact  that  the  solution  applies 
only  very  close  to  the  front  of  the  applied  load.  This  is  demonstrated  by  Figure  10, 
which  shews  also  the  surface  pressure  that  occurs  at  the  point  with  the  same 
abscissa  5  as  the  target  point,  and  at  the  time  of  arrival  of  the  plastic  front 
at  the  target.  This  surface  pressure  decreases  in  the  same  range  to  about  25$  of 
its  original  peak  value.  (It  is  noted  that  the  decrease  of  peak  pressure  with 
depth  for  the  decaying  pul3e  is  most  pronounced  for  small  values  of  ct  and  v,  and 
that  the  example  selected  shows  a  faster  decrease  of  the  peak  pressure  with  depth 
than  all  but  one  of  the  other  cases  listed  in  Tables  III  to  XI.  Tin  exception  is 


the  nonrepreoentative  case  shown  in  Table  XII.) 


Figure  10  also  shows  the  decrease  of  the  peak  pressure  with  depth  due  to 
the  decaying  pulse  for  an  elastic  material  having  the  same  value  of  Poisson's 
ratio  v.  There  is  a  slight  decrease  between  y  «  0  and  0y  ~  0.055,  but  for  larger 
values  of  y  there  is  no  further  decrease.  The  maximum  decrease  is  about  7$. 

The  effect  of  the  decay  of  the  pulse  for  the  elastic-plastic  material,  while 
quite  modest  is  several  times  larger  than  for  the  elastic  material. 

2.  Typical  results  for  V/cp  »  5 


Consider  again  the  case  v  »  0.125,  ®  -  0.10,  for  which  the  results  are  listed 
in  Table  III  and  the  accompanying  Fig.  11.  The  major  difference  between  the  results 
for  v/cp  a  2  and  those  for  V/cp  a  5  is  the  much  smaller  violation,  n  »  1.02,  near 
the  origin  for  the  latter,  and  the  fact  that  there  are  no  further  violations  at 
all  on  the  S-front  at  any  of  the  points  considered.  The  validity  is  limited  by 
the  violation  at  point  15  on  the  surface.  It  is  typical  for  this  value  of  V/cp 
that  violations  on  the  S-front  do  not  control,  but  that  violations  on  the  surface 
or  the  truncation  limit  the  validity  of  the  solution. 


It  is  useful  to  compare  the  results  for  V/cp  »  5  with  a  simpler  analysis  for 
the  limiting  ease  V/cp  -*■  »  .  To  consider  the  limiting  case,  the  definition 
§  =  x  -  Vt,  Eq.  (32 )>  is  introduced  into  the  surface  load  given  by  Eq.  (h6)  for 
x  <  Vt 


y 


-poe’ 


ue 


■  poe- 


a(x  -  Vt) 


lou; 


Consider  now  the  time  history  cf  the  surface  pressure  at  a  point,  say  x  =  0,  for 
V  -*  »  .  In  the  limit  ^  must  go  to  zero,  u  -*  0,  in  such  a  way  that  |*V  becomes  .the 
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inverse  of  the  decay  time  of  the  pressure, 

!*V  -  (81) 

o 

The  surface  pressure  in  the  limiting  case  is  then 

"  t/t 

y -p0«  (82) 

V  -*>  «•  implies  that  the  pressure  is  applied  simultaneously  at  all  points  on  the 
surface  of  the  half -space,  i.e.,  the  problem  becomes  a  simple  case  of  one 
dimensional  propagation  of  a  plane  wave.  The  case  is  treated  in  Appendix  III 
using  the  same  general  approach  as  in  Sections  II  and  III.  Using  this  formulation 
for  y/cp  h>  *»  ,  the  decrease  of  the  peak  pressure  with  depth  was  computed  for 
several  values  of  v  and  a.  The  results  for  the  peak  pressure  agree  to  two 
significant  figures  with  those  for  V/Cp  ■  5>  providing  a  check  on  the  two 
dimensional  analysis.  The  comparison  of  the  two  cases  also  proves  that  estimates 
for  the  magnitude  of  the  difference  between  the  effect  of  step  and  decaying 
pressures  for  values  of  V/cp  >  4  may  be  based  on  one  dimensional  analysis. 

3*  Range  of  depth  for  which  results  apply 

The  previous  discussions  considered  the  validity  of  the  results  only  with 
respect  to  a  nondlmensional  depth  j*y.  It  is  appropriate  to  consider  a  value  p  of 
physical  interest  to  find  out  for  what  actual  depth  y  results  have  been  obtained. 

hot  tQ  ■  0.2  sec.,  V  ■  4000  ft/sec.,  and  Bj.  (8l)  becomes 


i  ■  Vt  »  800  ft 
Ifr  o 


(83) 


Table  TI,  representing  a  typical  case  for  V/cp  *»  2,  contains  information  on 
the  peak  pressure  up  to  about  |&y  «  0.6  or  to  a  depth  of  about  t8o  ft,  A  time 
history  can  be  plotted  from  points  8  and  9  and  from  an  interpolated  point  on  the 
S-front  between  points  IP  and  13  for  ny  *  0.4  or  at  a  depth  of  300  ft. 

For  V/cp  a  5,  where  results  ore  r/hown  in  Table  III,  the  equivalent  results 
reach  only  a  depth  of  about  160  ft.  Howevei*,  in  this  case  the  one  dimensional 
analysis  in  Appendix  III,  valid  for  somewhat  larger  depth,  can  be  used  below  the 
depth  of  l6o  ft,  or  one  can  use  relatively  simple  finite  difference  schemes  to 
reach  any  depth. 


SECTION  V 


CONCLUSIONS . 

An  analytical  approach  has  been  presented  and  numerical  results  have  been 
obtained  for  the  effect  of  a  decaying  pressure  pulse  moving  with  superseismic 
velocity  V  >  2  Cp  on  the  surface  of  a  half-space  of  elastic-plastic  material 
subject  to  the  yield  condition,  (1) .  This  yield  condition  applies  to  a  material 
with  internal  friction  of  the  Coulomb  type.  Due  to  simplifications  in  the 
analysis  the  results  obtained  are  approximate,  and  restricted  to  target  points 
of  limited  depth,  as  discussed  at  length  in  Section  IV. 

Information  on  the  stress  field  for  exponential  decay  of  the  applied  load 
is  contained  in  Tables  2-12  and  in  the  accompanying  figs.  (10)-(20).  Comparing 
the  effects  of  step  loads  given  in  [l]  with  the  new  results  for  decaying  loads, 
it  is  found  that  in  the  latter  case  the  peak  and  subsequent  pressures  decrease 
with  depth,  which  is  not  the  case  for  a  step  load.  However,  the  decrease  found 
at  any  depth  y  is  only  a  fraction  of  the  decrease  of  the  applied  pressure  at  the 
point  oa  the  surface  directly  above  the  target  point,  i.e.,  a  distance  equal  to 
(-  y/tg  $)  behind  the  shock  front.  The  decrease  is  therefore  less  important  than 
one  might  eiqpect. 


It  has  also  been  found  that  for  velocities  V  >  4  c„  useful  information  for 

<v  r 

the  two  dimensional  ccse  of  progressing  surface  pressures  may  be  obtained, 
approximately,  from  one  dimensional  analysis,  which  is  naturally  much  simpler. 
E.g. ,  it  is  not  too  difficult  to  formulate  a  finite  difference  analysis  treating 


fho  n»*rthl  A  u  a  _ 

- - - - - — *|fl.  n^JA5UUi.A  V  'XLAVWXIJg  1  UX  MU?  ttUQi  1/1 01184. 

plastic  regions  ignored  in  the  Appendix.  The  validity  of  results  obtained  in 
such  a  fashion  would  th.jn  not  be  subject  to  the  depth  limitations  which  restrict 
the  results  obtained  here. 
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APPENDIX  X  -  Spooled  Case  of  the  Step  Load,  a  ■  0 


Although  an  approach  analogous  to  that  used  for  a  >  0  may  also  be  used 
for  the  special  case  Of  the  step  load>  o  ■  0,  it  is  simpler  to  take  the 
solution  found  in  Appendix  D  of  Ref.  [l]  directly.  As  discussed  in  Section  II 
the  solution  he.a  an  inelastic  pressure  front,  but  is  otherwise  entirely  elastic. 


#;'v4A| 


The  stresses  in  the  region  $  <  f  <  <pg  are 


.v  i 


where 


»  Aa 

o2  -  -  RAcr 


cos  apg 

A<r  -  »  p  - - g""V . . . 5 - 

(1  -  R)  cos  (<p  -  9S)  +  (1  +  R)  cos  9g  -  1 


1  -  Of/; 


and  o^  ,  Og  ,  o^  are  the  principal  stresses.  The  stresses  in  the  region 
tr  <  cp  <  9^  are 


=  -  P0  (88] 

®2  *  P0  +  (i  +  R)  Aa  (89] 

<»2  "  RAa  (90  J 

The  invariants  required  to  check  the  violation  of  the  yield  condition  are 

Jx  -  (1  +  2R)  Aa  (91) 


J2  m  Pq  +  (1  *  R)  poA~  + 


m  i 
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APPENDIX  H  »  Special  Case  of  the  Linear  Term,  n  »  1 


The  equations  governing  the  response  of  the  system  to  a  linear  surface 
loading  are  obtained  from  the  general  expressions  derived  in  Section  III  by 
setting  the  index  m  «  1.  The  potentials  are  then 


•  %y3  +  ftjSy2  +  tt2S2y  +  »3?3 
*{1)  -  b0y3  +  bjy2  +  bgS2y  + 

Y^  m  cQy^  +  cjy2  +  cg52y  +  c^53 


1 

(93) 


i 


There  are  twelve  unknown  coefficients  a^  ,  b^  ,  c^  ,  (i  ■  0  to  3),  and 
Eqs.  (71)- (79)  written  for  m  «  1  provide  twelve  equations  for  their  solution 


6(II  "  2)  a3  +  “  -  T1*  (9*0 

^2  "  "  2)  c3  "  0  (95) 

tan2«pg  (bx  -  cx)  +  2  tan  9S  (bg  -  cg)  +  30>3  -  c3)  »  0  {96) 

6Nja3  +  2(N1  tan  9  +  Kg)  ag  +  2Ng  tan  9  +  6N^b3  + 

+  2(N^  tan  9  +  N^)  bg  +  2N-  tan  9  b^L  »  0  (97) 


6N^a3  +  2(Nj  tan  <p  +  Ng)  8g  +  2Ng  tan  9^  +  6Ngb3  + 

+  2(Ng  tan  9  +  J^)  bg  +  2It^  tan  9  «  0  (98) 


6N^a^  +  2(N^  tan  $  +  &^)  ag  +  2K^  tan  $  a1  +  6N1Qb3  + 

+  2(N10  tan  9  -  Ng)  bg  -  2Ng  tan  9  bx  »  0  (99) 
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It 


.1 


1- 


* 


ax  »  30f  -  1)  a3 

a  .  &Ull 


*x  -  3(l|  -  1)  b, 


3  a2 


bo* 


(4  "  V 


3  ~  b2 


C1  -  3(^  -  1)  c 

<^-l) 


3  c2 


\ 


}  (100) 


Substituting  Eqs.  (100)  into  Eqs.  (  9k)  to  (99  )  reduces  the  system  to  six 
equations  on  the  six  unknowns  V  ^  .  V  bg  ,  Cg  .  ^  (  94) 

and  (  95  )  do  not  change,  while  Eqs.  (  96  )-(  99  )  become 


6(b3  “  c3)  +  2  tan  <*g  (b2  -  c2)  ®  0 

C6N1  +  6N2  taa  $  (M2  ”  1)1  *3  +  2Cnx  tan  9  +  NgJ  dg  + 

+  C6N4  +  6n3  tan  9  (l^j  -  l))  b3  +  2[N4  tan  9  +  bg  -  0 

C6n5  +  6n6  tan  9  (^  -  1))  e3  +  2[N?  tan  9  +  N63  &g  + 

+  C6Nq  +  61^  tan  9  (K?  -  1)]  b3  +  2f»8  tan  9  +  1^1  bg  «  0 

[6n9  +  6i^  tan  9  (M2  -  1)]  &3  +  2il^  tan  9  +  ^]  &2  + 

+  [6 K1q  -  6N6  tan  9  (^  -  1)3  b3  +  2[NiQ  tan  9  -  ^3  bg  -  0 


(101) 


(102) 


(103) 


(104) 


1 


i 


t 


! 


APPENDIX  HI-  The  Limit! as  Case  V/cp  -  o. 

As  a  llts:. ting  case  for  very  large  values  of  v/cp  one  obtains  the  one 
dimensional  problem,  where  the  surface  load  is  applied  uniformly  everywhere 
on  the  surface,  and  decays  exponentially  in  time  (see  Pig.  A-l). 


The  applicable  relations  between  stress  rates  and  strain  rates  are 


•  a  »  _ 

t  *  t  *  e  *  0 

x  s  xy 


*  Sv  _  U+v)(l-2v) 


5r  *  ElttJ  ay 


Integration  with  respect  to  time  gives 


O  m  O  m  <3  +  f_(y) 

Z  X  1-v  y  5ljr ' 


(105) 

(106) 


(107) 

(108) 


where  and  f ^  are  arbitrary  functions.  Finally,  the  equation  of  motion  is 

(1Q9) 


sr 


The  solution  of  these  differential  e  {uations  may  be  expressed  J/;  terns 


of  a  single  potential  #  and  of  one  arbitrary  function  g^(y  • 


*y  +  g^y) 


v  * 


yt 


.MUaat 

y  l-2v  yy 


(110) 

(111) 

(11?) 


f 
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where  the  potential  #  satisfies  the  wave  equation, 


*yy  *  c2  #tt 

cp 


and  the  functions  f*(y>  and  gjfr)  are  related 

V 


f.  y  yt  . . 

H..J/  dy 


The  boundary  conditions  are 


(113) 


(114) 


at  y  »  0 


a  »  «  p  e 
y  *o 


-  &et 


m«0 


at  y  =«  ct 


ary  ■  -  pcv 


q  at  l.-I.0  ,/ii  .  gf 

X  1  +  asr/3  r 


The  Taylor  series  expansions  for  f  has  the  form 
•  m±2 

• .  Y,w 

m»0 


#(“> .  (at)1  y' 

1*0 


i  .jm+2-i 


(115) 


(116) 

(117) 


(118) 


with  n  +  3  unknowns  . 

Substituting  this  series  into  Eqs.  (115)-(117)  in  the  manner  tiscuc-seu  m 
Section  III  yields  m  +  3  equations 


'P  (m  +  2  -  Dim.  +  1  -  i) 


i+2  z2 


(i  +  2)(i  +  1) 


/lin) 


a  » 

Za 


(-  l)m  l-2v 
ml  4(l-v) 


(120) 
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The  reuniting  values  of  ft,  permit  computation  of  the  stresses  a  and  a  from 

x  y  z 

Eqs.  (112)  and  (108),  respectively,  where  the  function  f^(y)  is  determined 
by  Eq.  (117). 
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